% this is the script for ploting sectioned particles(circles) in 2D

clc;
clear all;
close all;

%Read circle data file
f_id=fopen('2D-section.dat');
XYR=fscanf(f_id,'%f %f %f\n',[3 inf]);
f_return=fclose(f_id);

X=XYR(1,:);
Y=XYR(2,:);
R=XYR(3,:);
numCircles=size(R,2);
INDEX=[1:numCircles];

maxR=max(R);
minR=min(R);
NumCircle=size(R,2);

%Draw circles

hold
for i=1:NumCircle
    [x1,y1]=scircle1(X(i),Y(i),R(i));
    patch(x1,y1,rand,'EdgeColor','none');
end;

%colormap (gray);
axis equal;
    rand;

